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Abstract 

We derive the microcanonical partition function of the ideal relativistic quantum gas 
with fixed intrinsic angular momentum as an expansion over fixed multiplicities. We 
developed a group theoretical approach by generalizing known projection techniques 
to the Poincare group. Our calculation is carried out in a quantum field framework 
and applies to particles with any spin. It extends known results in literature in that 
it does not introduce any large volume approximation and it takes particle spin 
fully into account. We provide expressions of the microcanonical partition function 
at fixed multiplicities in the limiting classical case of large volumes and large angular 
momenta and in the grand-canonical ensemble. We also derive the microcanonical 
partition function of the ideal relativistic quantum gas with fixed parity. 
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1 Introduction 



In a previous work [1], we derived the microcanonical partition function of 
an ideal relativistic quantum gas of spinless bosons within a quantum field 
theoretical framework. Therein, we confined ourselves to the microcanonical 
ensemble with fixed energy and momentum. Here, we will obtain the partition 
function of the microcanonical ensemble of an ideal relativistic quantum gas 
with fixed energy-momentum and angular momentum. It should be pointed 
out from the very beginning that, for a quantum system, energy-momentum 
and angular momentum can be denned simultaneously only in the rest frame, 
i.e. if P = 0. In fact, more rigorously, one should rather speak of fixing the 
maximal set of commuting observables of the orthochronous Poincare group 
10(1,3)^, including energy-momentum and two observables formed out of the 
Pauli-Lubanski vector. However, for the sake of simplicity, we will retain the 



more colloquial expressions of "energy-momentum and angular momentum 
conservation" or "microcanonical ensemble with angular momentum conser- 
vation" being understood that we are working in the system rest frame, where 
the Pauli-Lubanski vector is proportional to the intrinsic angular momentum 
or spin of the system. 

The microcanonical partition function of an ideal relativistic quantum gas with 
fixed angular momentum, in a full quantum field theoretical framework, has 
never been derived for particles with spin. It has been obtained in an essentially 
non-relativistic quantum mechanical framework in ref. [2] and an equivalent 
calculation has been performed in ref. [3] by using a projection method, yet 
only for spinless particles. With the same limitations, an expression with an- 
gular momentum conservation, but without momentum conservation has been 
derived by Koba [4] . The classical limit has been worked out in ref. [5] . In this 
work we will get over these limitations and derive the microcanonical parti- 
tion function of the ideal relativistic quantum gas of particles with spin by 
using a group theoretical approach including the maximal set of observables 
pertaining to space-time symmetries, i.e. energy-momentum, spin, and parity, 
in a consistent way in a quantum field framework. Like in ref. [1], we will write 
it as an expansion over fixed multiplicities. 

The most general expression will be used to work out special cases, such as 
the grand-canonical (fixed temperature) and the classical limit. We will show 
that those, more familiar, expressions, can be recovered, thus supporting the 
correctness of our approach. 

Notation 

In this paper we adopt the natural units, with H = c = K = 1. Space-time 
linear transformations (translations, rotations, boosts) and SL(2,C) transfor- 
mations are written in serif font, e.g. R, L. Operators in Hilbert space will be 
denoted by a hat, e.g. R. The only exception to this rule are field operators, 
which will be written as a capital ^, while c-numbers field functions with a 
small if}. Also unit vectors will be denoted with an upper hat, but they will 
always be referred to explicitely in the text. 



2 The projection onto irreducible Poincare states 

In order to calculate the microcanonical partition function of an ideal relativis- 
tic quantum gas fixing the maximal set of observables pertaining to space-time 
symmetries, the use of group theory is both essential and enlightening. The 
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statement of the problem in this framework has been developed in refs. [6,7], 
generalizing known projection operator techniques used for internal, rather 
than space-time, symmetries, and we will outline it here. 

In ref. [1] we showed that the microcanonical partition function of a relativistic 
system can be written as: 

n = e s = 5> v |PilM , (i) 

where \hy) are a complete set of states for the localized problem (the quantum 
field in a compact region) and is a projector operator which selects the 
quantum states with the fixed values of the physical observables; S is, by 
definition, the entropy If only energy-momentum is involved, the relevant 
symmetry group is T(4) and the projector reads: 

P l = ( 5 4 (P-P) (2) 

Although 5 4 (P - P) is not a proper projector because P 2 = aP where a 
is a divergent const antPl we will maintain this naming even for these non- 
idempotent operators relaxing mathematical rigour. 

Similarly, for the maximal set of observables, the projection Pj should be made 
onto irreducible states i of the full symmetry group, namely the orthocronous 
Poincare group 10(1,3)^. These states can be built by diagonalizing operators 
or combinations of operators of the Poincare algebra. The first Casimir is the 
squared four- momentum P 2 with eigenvalue M 2 . The construction of physical 
states depends on whether M 2 > or M 2 = 0. We will consider only the 
former case, i.e. massive particles, and stick to the convention of ref. [8] for 
the definition of spin, that we will briefly outline. 

First, the Pauli-Lubanski vector is formed: 

W, = -(1/2) £ W >P- (3) 

vpa 

J being the generators of the Lorentz group, with the following properties: 
[W„P„}=0 (4) 

per 

w-P = o 

1 This is owing to the fact that normalized projectors onto irreducible representa- 
tions cannot be defined for non-compact groups, such as the space-time translation 
group T(4) 



3 



In view of the first commutation relation, if \p) is an eigenvector of P, so 
is The restriction of W to the eigenspace labelled by four-momentum 

P is defined as W(P). Since W(P) ■ P = 0, this four- vector operator can be 
decomposed onto three orthonormal spacelike four- vectors n\{P) , n 2 (P) , n 3 (P) 
which form a basis of the Minkowski space with the unit vector P = Pj \f[P 2 )'- 

w{p) = Y i W i {P)n i {P) (5) 
i=i 

It can be shown that the operators: 

Si(P) = Wi(P)/M (6) 

form an SU(2) algebra and are the suitable relativistic generalization of the 
spin angular momentum. The third component Ss(P) can be diagonalized 
along with S 2 = —W 2 jM 2 which is a Casimir of the full group 10(1,3)^, with 
corresponding eigenvalues A and J(J + 1). With a suitable choice of rij(P), 
i.e.: 

m(P) = [P]ei [P] = R 3 (y?)R 2 (£)L 3 (0 (7) 

ti being the unit vectors of spacial axes and [P] being a Lorentz transformation 
bringing the timelike vector Pq = (M, 0, 0, 0) into the four-momentum P with 
polar coordinates £, 9, tp; the eigenvalue A has the physical meaning of the 
component of intrinsic angular momentum in the rest frame along the direction 
of particle momentum P. Thus, with the choice (7), A is the helicity in the 
rest frame and J is, by definition, the spin of the particle. Since, from eqs. (6), 
(5) and (3) Si(Po) = Ji, the spin operators in the rest frame coincide with the 
generators of the rotation groups. 

The inclusion of parity adds another discrete quantum number II to the max- 
imal set of eigenvalues. Altogether, an irreducible state of the this group can 
be labelled by a four-momentum P, a spin J, its third component A and parity 
II and the projector onto such state Pj can then be written as: 

Pj = PpjAn (8) 



Generalizing the known integral expression for compact groups, one can for- 
mally write the projector (8) as an integral over the Poincare group: 

1 r 

Ppjau = - dimz/y d/i(g z ) U£(g z )g z (9) 
2=1, n 

where \x is the invariant group measure, z is the identity or space inversion l~l, 
g z G 10(1,3)^, U v (g z ) is the matrix of the unitary irreducible representation v 
the state % belongs to and g z is the unitary representation of g z in the Hilbert 
space. Like for the energy-momentum projector (2), the expression (9) is not, 
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strictly speaking, a projector because the group 10 (1,3) ''"is non-compact and 
PpjAn = °PpjAn with a divergent constant. However, we will keep calling 
(9) a projector relaxing mathematical rigor. Working in the rest frame of 
the system, with P = (M, 0), the matrix element U^(g z ) vanishes unless 
Lorentz transformations are pure rotations and this implies the reduction of 
the integration in (9) from IO(l,3) T to T(4) ® SU(2) <g> Z 2 [6] (replacing SO(3) 
with its universal covering group SU(2)). The general transformation of the 
orthochronous Poincare group g z may be then factorized as: 

& = T(x)ZA = T(a;)ZL fi (OR (10) 

where T(x) is a translation by the four-vector x, Z is either the identity or 
the space inversion and A = !_„(£) R is a general orthochronous Lorentz trans- 
formation written as the product of a boost of hyperbolic angle £ along the 
space-like axis n and a rotation R depending on three Euler angles. Thus 
eq. (9) becomes: 



Ppjxn = I E / d 4 * / dA ^(T(x)ZA)* J(x)Z A 

= \^J^j d 4 ^ / dA ^n^(A)*t( a )ZA (11) 

where z — if Z = I and z — 1 if Z = n. In the above equation, the invariant 
measure d 4 a; of the translation subgroup has been normalized with a coefficient 
l/(27r) 4 in order to yield a Dirac delta, as shown later. Furthermore, dA is the 
invariant normalized measure of the Lorentz group, which can be written as 
[9]: 

dA = dl_ n (0dR = sinh 2 £d£ — dR (12) 

Air 

dR being the invariant measure of SU(2) group normalized to 1, £ G [0, +oo) 
and are the angular coordinates of the unit vector n. 

In the rest frame of the system, where P = Po = (M, 0), Lorentz transfor- 
mations A comprising any non-trivial boost (i.e. with £ ^ 0) would make the 
matrix element U^(A)* vanish. Therefore A reduces to the rotation R and (11) 
to: 

Ppjxu = I E tAi / d A x (2 J + 1) J dR e 1 ^ IF £&(R)' f(x) Z R (13) 



z=i 



where we now use the ordinary symbol D to denote matrices of SU(2) unitary 
representations. The irreducible representations of this group are now labelled 
by eigenvalues of its generators J 2 and J 3 which coincide with the spin J and 
its third component A defined above by means of operators Si(P). 
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Taking into account that [Z, R] = 0, we can move the Z operator to the right 
of R and recast the above equation as: 



PpjAn 



(2vr) 4 



1 



/ 




2 



5\P - P)(2J + 1) j dR D^R- 1 ) R 



i + nn 

2 



(14) 



II being the parity of the system. The Eq. (14) is indeed the final general 
expression of the projector defining the proper microcanonical ensemble with 
P = (M, 0), where all conservation laws related to space-time symmetries 
are implemented. The nice feature of the above expression is the factoriza- 
tion of projections onto the energy-momentum P, spin-helicity J, A and parity 
II which allows us to calculate the contribution of angular momentum and 
energy-momentum conservation separately. 



3 The full microcanonical partition function 

The microcanonical partition function (MPF) (1) can be written also as, by 
inserting a resolution of identity J2f\f)(f\'- 



where Py = J2h v \hv)(hv\ is the projector onto localized states \hy) and 
Pi = PpjAn defined in Sect. 2. In ref. [1], we argued that Py can be written, in 
its most general form, in a field theoretical framework as a functional integral 
over eigenstates of the (scalar charged) quantum field: 



where the field function ip(x) is some arbitrary function outside the region 
In the functional integration (16), the index V implies that integrated 
field variables are just those in V, i.e. T>ip = Ilxgv dip{x). With this definition, 
spurious terms depending on the degrees of freedom outside V are introduced 
which must be subtracted "by hand" in a consistent way. This means that no 
contribution arising from these degrees of freedom should be retained at the 
end of the calculation [1]. 

2 We will use the same symbol V to denote the system region and its volume. 



^ = E(/l P ^ P ^l/)=tr[P i Py] 



(15) 




With |V, V f ) = ®xev|^(x),V t (x')> ®*{V ^(X^V)) 



(16) 
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The quantum field theoretical approach was suggested in ref. [7] to be essential 
for small volumes, because localized states are in principle a superposition of 
states with different number of particles (meant as asymptotic states of the free 
field defined over the whole space). The approximation of assuming localized 
states having a well-defined number of particles is expected to be a good one 
only when the size of the region is much larger than Compton wavelength. 
Yet, the final expression of the microcanonical partition function obtained in 
ref. [1] only differs by that calculated with the aforementioned approximation 
by an immaterial constant factor, after due subtraction of spurious terms. In 
this work, we will extend this result to particles with spin, that is to non-scalar 
fields. 

As a first step, the MPF Q is expressed as an expansion over all possible 
channels. By channel we mean a set of particle multiplicities for each species, 
that is a set of integers N 1 , . . . , N K = {Nj} where K is the total number of 
species. In formula: 

n = E %J • (17) 

Q{n •} is defined as the microcanonical channel weight and can be calculated in- 
tegrating on kinematical degrees of freedom. Writing |/) = |{A^}, {k}) where 
by {k} we label the set of kinematical variables (momenta and spins) of the 
free Fock space state |/), the microcanonical weight of the channel {Nj} reads: 

%v,} = E W TO.W = E<{^M*}l p i p vlTO,{*}> • (18) 
{k} {k} 

where u>{N },{k} is defined as the microcanonical state weight. It is now useful to 
insert a resolution of the identity in W{N-},{k} by using again the completeness 
of a set of Fock space states \{Nj}', {k'}): 

w = E ({Njh {kmmy, {k'})({N 3 y, {k'}\p v \{N 3 }, {k}) 
w >',{*'} 

(19) 

and work out the matrix elements (f\Pi\f) and (f'\Pv\f) separately. It is 
easy to realize that P; in (19) cannot change the multiplicities of the ket \f). 
This can be seen from its integral expression (14): translation, rotations and 
reflections cannot change the number of particles of the state vector they are 
acting upon. Therefore, in the sum in eq. (19) N' = N and: 

W = E<TO» ik}\Pi\{Nih m)({Nih m\Pv\{Nih {k}) (20) 
{k'} 

We will start our derivation of the MPF by considering the projector Pj = 
Ppjx without parity fixing, which will be addressed in Sect. 10. According to 



7 



eq. (14), the reduced projector Ppj\ can be written as: 

Ppjx = 5\P - P){2J + 1) JdR D^R- 1 ) R . (21) 

We will first work out the channel with one particle and calculate Q 1 . This 
will serve as a basis for the calculation of the general microcanonical channel 
weight in Sec. 5. 



4 Single particle channel 

Henceforth, we will stick to the notation of ref. [8]. Single particle states are 
written as | [p] , cr) and fulfill: 

P\[p],a)=p\[p],a) and S 3 (p)\[p],a) = a\[p],a) (22) 

with normalization: 

(b] )( x|[g],T> = 5 3 (p-q)cW. (23) 

With [p] we denote the SL(2,C) (universal covering group of SO(l,3)i) matrix 
corresponding to the Lorentz transformation transforming po = (m, 0) into 
p = (e, p). This is similar to the notation of eq. (7), with the difference that 
now the notation is extended to SL(2,C) and expressions like [p] 1 " are now 
meaningful. In the previous equations, square brackets in |[p],r) warns us 
that, strictly speaking, the corresponding state (and also the operator S^p)) 
has been defined relatively to a non-unique set of axes p, rij(p); i — 1,2,3 (see 
Sect. 2), but this will not affect our calculation in any way [8] and we will 
drop henceforth this notation letting: 

\\p],a) = \p,a) . (24) 

The transformation of a state |p, a) under a general Lorentz transformation A 
reads: 

N\p,o) = £ |Ap,r)Df CT ([Ap]^A[p]),/^ . (25) 

T V P 

where S is the spin of the particle, (Ap)°, p° are respectively the time com- 
ponent of Ap, p and [Ap] _1 A[p] is the Wigner rotation. A special case of (25) 
occurs when A is a rotation: 

R\p,a) = J2\^P,r}D s T(T ([Rp]- 1 R[p]) . (26) 

T 

According to eq. (20), we have to evaluate: 

= EMPibVXpVlPyb,^ . (27) 

p',cr' 
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By using the integral expansion of the projector Ppj\ in eq. (21), we can easily 
provide an integral expression of the first matrix element in eq. (27): 

(p,cr\P PJX \p > ,a > }=5\P-p)(2J+l) J dR Z^R" 1 ) 

xS 3 (Rp'-p)D s aal (Wr 1 m)- (28) 

Conversely, it is much more difficult to develop the expression of the second 
matrix element (p, <r|Py|p, a'). For spinless particles, this requires a functional 
method, described in detail in ref. [1]. To cope with this problem, we first need 
the expression of fields associated to particles with generic spin S. The free 
fields corresponding to the 2{2S + 1) physical degrees of freedom of a particle 
and antiparticle with spin S in Heisenberg representation can be written as 
[8] (see also ref. [10]): 



*r(x) = -irE /^^([p])«(p,^)e-^+^(MC'- 1 )6t(p, (T )e^ 
(27r)2 ; J V2e 

(29) 

where e = Vp 2 + m ' 2 is the energy; a, a) and b, are respectively destruc- 
tion and creation operators for particles and antiparticles; D s is the finite- 
dimensionsal representation of SL(2,C) of dimension 2S + 1 and the SL(2) 
matrix C = icr 2 (<7j being a Pauli matrix). The (anti) commutation relations 
of creation and destruction operators in (29) are related to the normalization 
of states (23) and read: 

[o(p, a), a\p\ a')] ± = [b(p, a), b\p>, a')} ± = 5 3 (p - p')5^ . (30) 

Under a Lorentz (or SL(2,C)) transformation, the fields in eq. (29) transform 
according to [8,10]: 

m T {x)K-^ =^Df ff (A)ti(Ax) . (31) 

a 

where we have written them in the Heisenberg representation. Note that D S (A) 
and D 5 (A)'I'~ 1 are two inequivalent finite-dimensional non-unitary representa- 
tion of the Lorentz group and correspond to those labelled as (S, 0) and (0, S) 
respectively. As representations of SL(2,C) they have useful properties: 

D S (A ] ) = D s (Ay D S (A T ) = D S (A) T . (32) 
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The matrix C = ia 2 fulfills following relations: 

C = -C T = C* C 2 = -1 CAC' 1 A T = (det A) I MA e SL(2, C) 
and plays a crucial role in the spin-statistics connection, as: 



(33) 



D»AC 2 ) = (-lf S 5, 



(34) 



The 2(2S + 1) field degrees of freedom in eq. (29) are needed to represent par- 
ticle and antiparticle. Since *& T (x) = D^ a (C)^^ (x); where \l/ c is the charge- 
conjugated field which is obtained from \1/ by swapping a <-> b and a* <-> 
6' in the eq. (29), for a particle coinciding with its antiparticle, \I/ T (x) = 
D^(C)^(x), a constraint which effectively halves the field degrees of free- 
dom. 

It is easily seen from (29) that for equal times t — t' [^(x), ^(x')] ± = 0, where 
the sign + or — applies to the fermion or boson case respectively, so it is 
possible to put together the fields \1/ and \I/ in one spinor with 2(2S + 1) 
independent (anti) commuting components: 



v|/ = 



(35) 



and write the two fields in eq. (29) in a compact way as: 
1 ^ r d 3 p 



v|/(x) = 



E 



C/(p, cr)a(p, a) e"*-* +V(p, a)^, a) e ip x (36) 



(2tt)3 V2e 
where U(p, a) and a) are 2(25' + 1) -dimensional spinors: 



C/(p,<7) 



^L([p] 



V(p,<7) 



(37) 



It is also possible to write the field in an even more compact way defining the 
2(2 S + 1) x (2S + 1) matrices U(p) and V(p) by clumping the spinors U(p,cr) 
and V(p, a) respectively: 



u(p) = ( DS([P]) 



V(p) 



(38) 
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so that the eq. (36) can be recast as: 

^ = 1 4r ( u (p)°(p) e ' ip ' x +v(pP(p) (39) 

(27r)2 J y/2e v 7 

ct(p) and b'(p) being 2S + 1-dimensional column vectors of destruction and 
creation operators. Now, let us define the (2S + 1) x 2(25* + 1) conjugate 
spinorial matrices: 

U(p) = U\p)T V(p) = V\p)F (40) 
r being the 2(25 + 1) x 2(25 + 1) matrix: 




The following relations apply: 



U(p)U(p) = (-l) 2S V(p)V(p) = 21 



U(p)V(p) = V(p)U(p) 



j2D s (C) if S even 
[0 ifSodd 



(41) 



as well as: 



U(p)U(p) 
V(p)V(p) 



I 



I 



-1) 2S I D s ([p}[p]^ 



dHIpMT 1 (-i) 25 i J 

One can construct the conjugate field W: 

d 3 p 



(42) 



(27T) 



2e 



(a f (p)C/(p) e ip - x +6(p)F(p) e"*-*) (43) 



The fields v|/ and W obey homogenous generalized Dirac equations for spin S 
[8,10]. Particularly, for S — 1/2 it can be shown that ^(m)U' is the usual Dirac 
field and v|/ obeys the familiar Dirac equations in the Weyl representation: 



i^d^ - = 7° = r Y 



-a 1 



(44) 
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To proceed with the calculation of the matrix element (p, a\Pv\p, &') we now 
have to express the creation and destruction operators as a function of the 
fields. We do this by using the fields in Schrodinger representation, i.e. with 
the fields at t — 0: 

(0|o(p,(7) = (0| — / d3x e-*P-v^it7(p) OT V T (x) 
2(27r)2 t j 

gW)|0> = -^-t£ / d 3 x e^v^M/^x)^)^ |0) 
2(27r)2 t j 

(0|6(p,(7) = (0| -^-tE / d 3 x e-^V2e^ T ^)V{p) TCJ {-lf s 
& t (P^)|0) = -^-rE / d3x ^^V(p)A(x)(-l) M |0) ; (45) 

these relations can be checked on the basis of (41). Therefore, using (45) for 
particles: 

(p',a'\P v \p,a) = (0|a(pV)PyaW)|0> 

= |^ / dV/ d 3 x e-^- x '^ x Ef/(p')^(0|^(x')Pv^(x)|0)f/(p) rCT 

(46) 

The difficult task in the above equation is to calculate the vacuum expectation 
value in presence of the projector Py. For the scalar case, we could cope with 
this problem through a functional integration of the field degrees of freedom 
in the region V and dropping the terms involving the points outside the region 
V [1]. However, the functional method developed for the scalar field cannot be 
trivially carried over to the non-scalar case for a twofold reason: firstly, even 
in the bosonic case, the fields W and vj/ do not commute at equal times and, 
as a consequence, writing an expansion of the operator Py like in eq. (16) is 
not possible; secondly, one has to deal with fermionic degrees of freedom. 

Yet we have some clues of how to work out (46), getting inspiration from the 
scalar case. The presence of Py means that eventually we have to restrict all 
spacial integrations to the region V. Then, we will assume that: 

/ d 3 x' e -P'- x '(0|vMx')Py^(x)|0> = — ([/(p')[?(p')) r , r <0|P y |0> (47) 
and: 

/ dV e-P'- x \0|vj/ T ,(xOPyV|/ T (x)|0) = ^^(y( P OnpO)r'r(0|Py|0). (48) 
Jv le' 

The justification of the formulae (47) and (48) resides in two facts: 
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• they are consistent with the limit Py — > I, as it can be readily checked from 
the explicit field expression and the normalization condition (23); 

• they have been proved for the scalar case [1]. 

If one accepts the conjecture expressed by the equations (47) and (48) the 
eq. (46) turns into, by using (41): 

(p',a'\P v \p,a) 

= l^F v (p-p>)(U(p>)U(p)) a , a (0\P v \0) 
2 V e 

= \ \[^MP - p') (^([PTW + V" 1 ))^ (0|Pv|0> (49) 

where F v is a Fourier integral over the system region V: 

F V (P ~ P') = j^y 3 J y d 3 x e--(P-P') . (50) 

For p = p' and a = a' the above expression reduces to V r /(27r) 3 (0|Py|0) which 
is the same we obtained for the scalar case in ref. [1]. It is straightforward to 
show that the same expression holds for antiparticles. 

We are now in a position to calculate the microcanonical state weight uj\. By 
plugging (28) and (49) in the definition (27): 

= 5\P- P )j:jd 3 p' (2J+l)jdRD( x (R' 1 )5%R P '-p)Dl l ([Rpr 1 R[p']) 

a 1 

x 1\^mp - p) (d'upvip]) + jaipiw- 1 )^ <oipv|o> 

= 5 4 (P-p)£/ d V (2J+l)JdRDi x (R^)5%Rp'-p)Dl l ([p]- 1 R[p']) 

a' 

*\Mp - p') (D^ip'rip]) + D^lp'W 1 )^ (0|Py|0> (51) 

Two matrix products are to be computed in the above equation. The first 
yields: 

E^(W" 1r Ip1)^(M" 1 W)) = DS (\pr lR \p]U (52) 

a' 

To work out the second matrix product we observe that [Rp'] _1 R[p / ] = [p] _1 R[p'] 
is a Wigner rotation, hence it is unitary: 

[p]- i RW] = ([prwir 11 = (IpVr-w = (Mwr 1 ) 
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thus: 



E^(M" 1 R^)^([pTM t ' 1 ) = E^([p] t RM t ~ 1 )^(bTb] t " 1 ) 

cr' a' 

= Dl(lpMpt 1 ) (53) 
By using (52) and (53) the eq. (51) becomes: 

Wl = 5\P -p)J d 3 p' (2 J + 1) / dR D^R-^Rp' - p) 

x ±F v (p - p') (^([pJ^Rb]) + ^([^R^- 1 ))^ (0|P y |0) (54) 

This is the final expression of the single-particle channel. We note in passing 
that, with the same method, the more general matrix element (p, (r'\PiPy\p, cr) 
could have been calculated, yielding: 

(p,a'\P t P v \p,a) = 5 4 (P-p) |d 3 p' (2J+1) /dRZ^R-^W-p) 
x±F v (p - p') (D^lp^Rlp}) + ^(^tR^t-i))^ ( |P y |0) (55) 



5 Multiparticle channels 



We first tackle the problem of a channel with N particles of different species. 
Writing the Fock state \N, {k}) in eq. (20) as: 

\f) = \N,{k}) = \N, K},K}> 

where {p n }', Wn} label the set of four-momenta and helicities of the N parti- 
cles, the microcanonical state weight in eq. (20) becomes: 



N "I 

II / d 3 P ; (N, {p n }, {a n }\P PJX \N, {p'J, {a' n }) 

.71=1 J 



{<} 

x(N, {p' n },{<}\Pv\N, {Pn},K»- 



(56) 



Both matrix elements of Ppj\ and Py can be calculated by a straightforward 
generalization of the single-particle channel because, as particles are distinct, 
creation and destruction operators, as well as field operators, commute with 
each other and the problem is fully factorizable. Therefore: 
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(N, {p n }, {a n }\P PJX \N, {p'J, {<}) = <5 4 (p - f; Pr )j 



x (2J+l)|dR^ A (R- 1 ) 



N 

n 

Ln=l 



n^dw'Ripjj^Rp^-pn) 



(57) 



and: 



W \Pn}, Wn}\Pv\N, {p'J, {<j n }) = (0|P V |0> 



X 



JV i f— - 

II oXhrMPn - p'J (D'UriPn]) + DSWW- 1 ) , ■ (58) 

Putting (57) and (58) into (56) and working out the matrix products, like in 
eqs. (52), (53) we get: 



r r N r 

^=(0|P y |0)(2J+l) /dRD^R- 1 ) n /d 3 pU 3 (Rp;-Pn) 

J Ln=l J 

1 -i / W 

x - (^([pJ^RH) + ^([pJtR^Jt- 1 )) F y ( Pn - p'J <5 4 P - £ p. 



(59) 



The integration over momenta p^ can now be performed. Since p' n = R p n : 

W/ = (0|Py|0) 5* (p - f>„) (2 J + 1) / dR D^R" 1 ) 
AT i 

II o (^ 5 ([P«]" lR W) + ^([z^Rk]*- 1 )) Fy(Pn - R _1 p 



(60) 



which is the final expression of the microcanonical state weight for N different 
particles. Integrating (60) over momenta p n and summing over spin projections 
<r n , one gets the microcanonical channel weight in a simple form: 



Q N = (2J + 1) /dR -Daa(R _1 ) II / d V trD 5 "(R) 
x F v {p n - R-'pJ 5 4 (p - J2Pn) (0|Pv|0) . 



(61) 



n=l 



where we took advantage of cyclicity property of the trace. It is interesting to 
note that Wigner rotations and momentum-dependent matrices in the matrix 
element of Py (58) combine so as to leave, once the sum over spin projections 
is carried out, a trace of a simple SU(2) rotation. 
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We will now deal with the more complicated case of identical particles. 
5. 1 Identical particles 

Developing the matrix element (/'|Pv|/) f° r a multiparticle state with identi- 
cal particles is far more difficult than for different particle species because the 
associated fields do not commute and factorization does not hold; this has been 
shown for the scalar field in our previous work [1]. Nevertheless, we showed 
that, keeping only terms involving field degrees of freedom in the region V, 
the result is the same one would get in a multi-particle non-relativistic quan- 
tum mechanical framework enforcing symmetrization or antisymmetrization 
for bosons or fermions respectively. Although the conclusion was fair, the proof 
involved lenghty calculations which would become even lenghtier for particles 
with spin. Therefore, we will not tackle the proof here and simply assume that 
it holds for particles with spin: a general multiparticle state can be treated as a 
state with different particles, with (anti) symmetrization for identical particles. 
Accepting this conjecture allows us to work on the fixed-number multiparticle 
tensor space instead of the original Fock space. 

Let pj be a permutation of the integers 1, . . . , Nj, x(Pj) & s parity and bj = 0, 1 
for bosons and fermions respectively. The generic state |/) in the multiparticle 
tensor space can be written as: 

ni^'K-K)},KK)})- (62) 

where K is the total number of particle species and {p n .}, {a n } the set of 
momenta and spin projections of particles of species j in |/). Also let p be the 
set (pi, . . . , pk) of permutations of integers (1, . . . , Ni), . . . , (1, . . . , Nk) and 
let N = J2j=i Nj. Let us also introduce a shorthand \f p ) for the ket: 

l/p> = ni%{Pp>,)M^K)}> (63) 

which is the multiparticle tensor formed by exchanging particle indices by the 
permutations p. By using this notation, the state |/) in (62) can be rewritten 
as: 

l/HE 
p 

and the microcanonical state weight uif. 

</„|PiPv|/p> (65) 
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yr x(Pj 

3=1 



x(pj) 1 

3=1 



l/p 



(64) 



w / = EE 



K 

n 

3=1 



xiPjVj] 
Nj\ 



which can be simplified, redefining dummy particle indices, into: 



n x(p 3 y 



(MPiPvlfp) 



(66) 



where t is the identical permutation. 



We first consider, for sake of simplicity, a multiparticle state with N identical 
particles. Eq. (66) reads, by taking = Ppj\: 



^f = Y,x(p) b (f.\Ppjx Pv\f P ) 

p 

= J2x(p) b (N, { Pn }, {a n }\P PJX P V \N, {p p{n) }, {a p(n) }) 



(67) 



Like for eq. (19), we insert a resolution of the identity between Pj and Py 
leading to: 



LIE/ dX 

1=1 cL 



(N, { Pn },Wn}\Ppjx\N, K},K}) 



x (N, {p'J, W n }\Pv\N, {p p(n) }, {a p{n) }) 



(68) 



Both matrix elements of Py and Ppj\ can be evaluated, in the multiparticle 
tensor space, the same way as for distinct particles. By using eqs. (57) and 
(58): 



W/ =(0|Py|0>5 4 (p-f>n) ITE/ d3 Pn 
\ n =l / n=l a' J 



(2J+1) 



U^Pn-P^DlAiRpX^Wn}) 



p J 

X ^MPp(n) ~ , 

V E n 4 °Vi°>(n) 



(69) 



If we now carry out the integration over momenta p\ and sum over all spin 
projections a' n we get: 



17 



"f = S 4 (p-j: Pn) (2 J + 1) / dR ^(R- 1 ) £ X (P) 

\ n=l / J p 

- N -j 

II 5 (^(M^RbpwD + ^dPnl+R^n)]*- 1 

|Pv|0>. 



.71=1 



l £ p(n) 



Fv(Vp{n) ~ R 1 Pn) 



(70) 



Indeed, the factor Il^Li \j £ p{n)/£n can be dropped from eq. (70) because nu- 
merator and denominator are altogether equal being p a permutation. 

The microcanonical channel weight of iV identical particles can be calculated 
by integrating Uf in (70): 



^ = (oiPHo)( 2 j+i)/dR/j 3 : A (R- 1 ) £ n/ d3 p 

J ai,...,a N p iv - ln=l J 



b [ N 



X 



- (-D s ([p„] _1 R[p p ( n )]) + D 5 (b„]t R[ ^ (n) ]t-i)) Fy(pp(n) _ R -i pn ) 



X 



* 4 Kf>") 



(71) 



where the factor 1/N\ has been introduced in order avoid multiple counting 
of momenta and spin projections. 

The above formula can be further developed in the limit of large volumes, 
when: 

In this case, the two matrices of the irreducible representation of SL(2,C) in 
the formula (71) become equal. To show this, one first notes that the Dirac's 
delta implies p p ( n ) = R -1 Pn, hence: 



[Pn] -1 R[ZV(n)] -> k] _1 R[R _1 Pn] 

IPnl^tPpWl*" 1 - [Pnl^fR-V]^ 1 (72) 

The right hand side of the first equation is manifestly a Wigner rotation, 
hence it is unitary and equal to the right hand side of the second equation. 
This implies that in the large volume limit: 

\ (D s (\p n ]-iR\p p{n) }) + D s ([p n M Pp(n) ]^)) ^ D s ([ Pn ]^R[ PpH ]) (73) 

We observe that the same equation applies in the non-relativistic limit, where 
the matrix [p] is unitary, since the associated Lorentz boost in eq. (7) can be 
approximated with the identity. 
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An useful development of the formula in the thermodynamic limit can be 
achieved by decomposing permutations into the product of irreducible cyclic 
permutations, that is writing: 

p = c x ...c H . (74) 

Instead of working out the most general case, we will focus on an example 
which will show more clearly how to take advantage of the above decompo- 
sition; the generalization will then be straightforward. Let us then consider a 
special case, with 5 particles and the permutation: 

(1,2,3,4,5) — (3,1,2,5,4) (75) 

which can be decomposed in two cyclic permutations: 

ci = (1,2, 3) —►(3,1,2) and c 2 = (4, 5) — > (5, 4) . (76) 

If we now write down the corresponding matrices product in eq. (71), by using 
the reduction (73) we have: 



n=l 

x^ 1CT3 ([pi]~ 1 R[P3])^ 2 (N- 1 R[P2])^ 2(71 (N- 1 RN) 

where the first two factors on the right hand side correspond to the cyclic per- 
mutation C2, while the other three to c%. Summing then over spin projections 
a n , we end up with the simple and nice result: 



E A Di^JWr^lPpin)}) =trI> 5 (|p4]- 1 R 2 |p4])trI> 5 (|p 1 ]- 1 R 3 |p 1 ]) 

{ct„} n=l 

= trD 5 (R 2 )trZ> s (R 3 ) (77) 

In general, for each permutation p and its cyclic decomposition, if we let I be 
the number of integers in each cyclic permutation and hi the number of cyclic 
permutations with / elements in p such that E^i Ihi = iV0, one has: 

N N h ( \ 

e n d^j^rm) = n M 5 ( Ri )] ( 7§ ) 

{a} n=l 1=1 

Thus, by using eq. (73) and its consequence eq. (78) in the large volume limit 
the microcanonical channel weight (71) turns into: 



3 The set of integers hi, . . . ,/tjy = {hi}, is usually defined as a partition of the 
integer N in the multiplicity representation. 
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Q N = (2J+1) JdR D a J a(R^)E 



x(p) 1 



xF y ( Pp(n) -R- 1 p„) [tr£> s (R") 



h n (p) 



n / d v 

n=l J 

S 4 (p-EPn) (0|P V |0>. (79) 



n=l 



an expression which also applies to the non-relativistic limit. On the other 
hand, if V is not large and we have a full relativistic regime, the two D s 
matrices in eq. (71) are no longer equal and the summation over polarizations 
gives rise to a very involved expression. The general (71) and the special case 
(79) can be easily extended to a multi-species ideal gas. They read: 



= (2J+l)/dR^ A (R- 1 ) £ £ n^li /d 3 p n 
J o-i,...,o-jv p |_j=i iy r n j =l J 



N 



X Fy( Pp{nj) - R-^.) 5UP-Y,Pn) (0|P V |0) 



n=l 



and: 



(80) 



( 2 j+i)/dR^ A (R- i )E n^n /d 3 P 

P 7=1 i' n ? =l 



hnApj) 



x F y (p, K .) - R" 1 ^) tr£> s '(R"') njW ' 5 4 [P- Y,Pn (0|Py|0) (81 



A' 



n=l 



respectively. In the eq. (81) Ylnj=i n jhj(Pj) = Nj. 

Finally, we observe that the eq. (81) could be further developed in the large 
volume limit so as to obtain a microcanonical channel weight which is re- 
summable analytically yielding a compact integral expression of the MPF, 
just like in the non-relativistic case without angular momentum conservation 
[6]. 



In the large volume limit expression (81) there appear traces of rotations. 
Since the trace is the invariant over group conjugacy classes, it turns out that 
all SU(2) transformations of the same angle ip around an axis n (which labels 
the different members of the same conjugacy class) have the same trace. It is 
therefore convenient to use the axis-angle parametrization for the integration 
measure over the SU(2) group: 
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fdR= — ^ r d9 d(f) dtp 2 sin 9 sin 2 $ (82) 
J lo7r^ Jo Jo Jo 2 

=i^/ dn *r d * 2si 4 

where (#, 0) are the polar and azimuthal angles defining the axis n. With this 
parametrization, the trace of a rotation Rn(V0 reads [9]: 



sin 

trD s (R h W)= \ - ^ ,J . (83) 
smf 



3in[(5+i)V] 
Furthermore, since: 

Rn(^l)Rn(^ 2 ) = Rn(^l +^ 2 ) => R&W = Rfi(^) 

the trace in eq. (81) simply becomes: 



sin [(5 + \)n^\ 



6 Spherical region 



The eq. (80) is the most general microcanonical weight of a channel {Nj} 
with the enforcement of energy-momentum and angular momentum conser- 
vation; eq. (81) is an approximate expression of (80) which applies for very 
large volumes. Summing (80) or (81) over all possible channels yields the full 
microcanonical partition function. 

If the system region has spherical symmetry, those expressions can be further 
simplified because the whole SU(2) group integrand in all the equations is 
independent of the axis n (see end of Sect. 5). The proof of this statement 
proceeds in two steps: first, one observes that flnv-} as a whole does not depend 
on the polarization index A for a spherical region, so that one can rewrite it 
summing over A and dividing by 2 J + 1. In the second step, one shows that 
the overall SU(2) integrand does not depend on the axis n. In formulae: 

n {N . } = (2 J + 1) / dR J D J (R- 1 ) AA /(R) = £ / dR ^ J (R^)aa/(R) 

A 

= J dRtr7J J (R- 1 )/(R) (85) 
where neither trD J (R _1 ) nor I(R) depend on the axis polar coordinates. 
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Let us prove the first part of our statement by writing the projector Ppj\ for 
a system at rest (i.e. with P = 0) as: 

P PJA = \M,J,\)(M,J,X\ (86) 

(where M is the mass of the system) and the channel weight in (18) as: 



^{N,} = E<{^}> J > A >< M > J > A|P y |{7V.}, {k}) (87) 

{k} 

= (M, J, X\P V (]T \{N 3 }, {kyUiNj}, {k}\) \M, J, A) . 

\w / 

Evidently, for spherically symmetric systems the condition: 

[Pv,R] = (88) 

must hold, because Py (see e.g. eq. (16)) involves the sum of all possible field 
configurations within a spherically symmetric region and the representation 
of a rotation on the Hilbert space transforms Py into itself. Likewise, the 
operator: 

Pw} = El{^}>M)({^}>MI (89) 

M 

is also invariant under rotations since it is a sum over all possible kinematical 
configurations of free states. Thus: 

[P{N J} , R] = (90) 

and, because of Wigner-Eckart theorem: 

n {Nj} = (M,J,X\P v P {Nj} \M,J,\) = (M,J\\P v P {Nj} \\M,J) , (91) 

as PyP{jv-} commutes with rotation operators according to (88) and (90). 
Therefore, the channel weight is independent of A. One can then sum over 
spin projections A the matrix element in all expressions of the microcanonical 
channel weight, and divide by (2 J + 1), i.e.: 

1 «r-^ r 1 r 1 srn ( J + \ ) ip 

pTTljp&W'M) - pjTij^WW) - (2J + 1) { ~f 

(92) 

where we used tr[R fi 1 (-?/')] = tr[Rfi(— ip)] = tr[R^('0)]. Thereby, the dependence 
on n owing to the matrix D J (Rn(ip)) is removed. 

The second step of the proof is to show that the function /(R) in the eq. (85) 
is also independent of the axis polar coordinates. By comparing eq. (85) with 
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the integral expression of Pja, i.e. 

P JA = (2J+1) J &RD J (R~ l )\\R 
and comparing (87) with (85) one has: 

W) = E({ N jh {k}\S 4 (P ~ P)RPv|{^-}, {k}) = tr[5 4 (P - P)RPyP {Nj} ] 
{k} 

(93) 

using definition (89). Let us now replace in the above equation the rotation 
R with axis n with that around the rotated axis On, which is ORO -1 . The 
rigthmost term in eq. (93) turns into: 



tr[5 4 (P - P)bR6- 1 P v P {Nj} ] = tr[05 4 (P - P)RP v P {Nj} b- 1 } 
= tT[5 4 (P-P)RP v P {N]} ] 



(94) 



where we used eqs. (88), (90) and took advantage of the commutation between 
S 4 (P — P) and rotations if P = 0. The eq. (94) proves our statement, namely 
the function J(R) is independent of the rotation axis and so is the whole group 
integrand function in eq. (85). 

In conclusion, for a spherically symmetric region, one can integrate away the 
solid angle fl^ and choose an arbitrary rotation axis, e.g. k, in all expressions 
of the microcanonical channel weight f2{jv-}- The integration measure of the 
SU(2) group (82) in the axis-angle parametrization effectively reduces to: 



2sin' 



(95) 



Particularly, the two equations (80) and (81) for a spherical region become: 



1 f An tb ( 1\ 

^} = — / d^ sin-sin J + -U ]T ]T 

Z7r J ° Z V Z/ <n,...,<r N p 



K X W n 

j=l ly 3- n 3 =l J 



n 



x \ (^(kJ^RsW^K.)]) + ^(W t RaWIPp(n J )] t - 1 



x F v ( Pp(nj) - R 3 (^)- 1 p„ J )l 5 A (p-jtpn) (0|Py|0) (96) 



n=l 



and: 
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n 



{Nj} 



2tt Jo 



r */> si,, i Bin n%£n/ ^ 

JU Z V Z/ p 7 = 1 iV i- n,- = l J 



sin(5j + 1)^-0 


h nj (pj)~ 


[ sin(^) j 





N 



^| P -EPn)(0|Pv|0) 

(97) 



n=l 



7 Partial wave expansion 



So far, we have been using a plane wave expansion to calculate the microcanon- 
ical channel weight. One may wonder whether an equivalent formula could be 
obtained using single particle states with definite angular momentum. In this 
section, we will show that this is the case. 

Instead of starting from scratch pur calculation by expressing the microcanon- 
ical state weight (19) with kets \jmSa) or \ jmlS) (in the /S'-coupling), we will 
work out the formula eq. (80) expanding plane waves into partial waves. For 
the sake of simplicity, we restrict to the simple case of Boltzmann statistics, 
corresponding to retain only the permutation identity in the sum in (80), 
yielding: 



n {Nj} = (2J+l) j dR D^R' 1 



n w n / 

7=i iV r m=i J 



x 1 [ A \V x ' ,p 

x (2tt)3 y 

= (2J+1) fdRD&R- 1 ) 



N 



n=l 



^tr£>*(R)j 5 4 ^P-J2Pn) (0|Py|0) 

N 

n / ^. 



x 



— *— / d 3 xe 4X -( p "- rlp "' trL> 5 "(R)l 5 4 (p-Y^Pn) (0|Py|0) 
( 27T ) J J V n=l J 



(98) 



where we used the explicit form of integrals (50) and performed the sum over 
polarization states in eq. (80) for the simple case of only identical permutation. 
In the second equality we dropped different indices for different particle species 
and used, for simplicity, a common index n for particles in the channels. 

The exponentials in Fourier integrals in eq. (98) can now be expanded into 
partial waves according to the well known formulae: 
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OO In 

e* x - p "=E E i'"47r^(p n x)^ mn (x)^ mn (p n ) (99) 

it 

OO * n 

e- lx ' Rp "=E E (-i)'"47r^(p n x)^ m ,(x)y,, m ,(Rp n ) 

where ji are spherical Bessel functions and Y are spherical harmonics; x and 
p are the unit vectors of x and p respectively. Consequently, Fourier integrals 
can be rewritten as: 



JL / d 3 X e^P»- R P") = (47T) 2 E ^"(-0^^„(Pn)^(RPn) 

v / l n ,m„; l' n ,m'n 

X (27p X d3x ^(P« X )x(Pn X )^rn„(x)F z r m , (X) . (100) 

Switching to spherical coordinates d 3 x — > x 2 dxdf2x, the angular integration 
can be made at once yielding: 

J d^ m (x)r; m ,(x) = s w 8 mm . , (ioi) 

hence eq. (100) turns into: 



>7r) 3 X d 



(27T) 

= j^p E ^(p»)^(Rp») X dx x2j ' (p " x) • (102) 



3 X e «x-(p„-Rp„) 



In eq. (98), for each particle, a factor trD 5 (R) appears. By taking advantage 
of the relation: 

Yi m (Rp) = E^'(P) D l m , m {RY (103) 

m' 

the product between the trace and Yi m (Rp n ) can be worked out as: 



tr/^(R)E>Wpn) £fc<(R) 

= E n.<(Pn) ^„(R)^<(R) 

= E ^nm^P") (^"' m n\ R^n; °Vii ^n; m n) 

m{„er„ 

= E ^^(Pn) (o , n,^|5' n Z n |J n ,^ n )(J n ,^|S , n Z n |<T n ,77l n )D^; Al{i (R) (104) 
m' n ,a n 
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where (mi, m,2\ji, j2\j, m) is a shorthand for the Clebsch-Gordan coefficient 
(ji, mi,j2, m 2 \j, m, j\,j2). In the last step, the rotation R has been expanded 
over a complete set of vectors according to: 

R = E E |J,/i>(^|R|JV)(JVI = E E l-MUVI i>i,(R)*jj' • 

(105) 

For each particle in the channel, a factor like (104) appears. The SU(2) group 
integral in eq. (98) can now be calculated collecting all factors depending on 
R: 



(2J+1) J dR ^-^^(^...^(R) 

= (2J + l)(J 1 ,n 1 ;...; J N ,n N \ j dR D^R- 1 ) R| Ji, /4; • • • 5 Jn, /4) 
— {Jii A*i j • • • '■> Jn-, Hn\P j\\Ji, Hi] ■ ■ ■', Jn, A* at) 

= (Ji, a*i; . . . ; J n ,^n\J, A)(J, A| Ji,//i; . . . ; Jn,v'n) ( 106 ) 

Collecting previous results, eqs .(99)-(106), the microcanonical channel weight 
(98) can be rewritten as: 



{Nj} 



3=1 3 



r N r 
n / d 3 P, 

_n=l 



N 



n=l 



X 



E (^i> A*i! • • • ; J/v, A*at| J, A)(J, A| Ji, . . . ; J N , fj,' N ) 



l n ,m' n ,<j n ,m n 
N 



X 



II - / dx « X ni/l(Pn X n) r /*mJPn)^nm;(Pn)K,^U^n/n|Jn,Ain) 
„=1 ^ J ^ 

x(J n ,/i^|S' n / n |o- n ,m n ) (107) 

This expression finally gives a clear physical meaning to our manipulations: the 
microcanonical channel weight of the channel is obtained by summing over all 
possible single-particle states \p n ,l n ,m n , S n ,<j n ), l n being the orbital angular 
momentum, S n the spin angular momentum, m n and a n their respective third 
components; and projecting these states first on states with total particle 
angular momentum J n and then all of them onto a state with total angular 
momentum J and third component A. 

A considerable simplification occurs if, in eq. (107) momentum conservation 
constraint is dropped. In this case the <5 4 reduces to 5 (m — J2n=i £ nj and 
the integration over the angular coordinates of three-momenta can be done at 
once yielding: 

/ d%^: m Jpn)*WP") = <W< • ( 108 ) 
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Thereby, the eq. (107) reduces to, using the unitarity of Clebsch-Gordan co- 
efficients: 



n 



{Nj} 



k I 

3=1 3 



_n=l J J V >2=1 



x E (^i,A«i;---;^Ar,A«Ar|J,A)(J,A|Ji,/x / 1 ;...; J N ,n' N ) (109) 
n 2 



X 



II - / y dX « X n jf„(PnXn) . 



Let us now restore momentum conservation. In the rest frame of the system 
where: 

( P ~ E Pn) = 6 ( M - £ e n ) f£ Pn) 

\ n=\ ) \ n=l ) \n=l / 

the <5 3 is expand into spherical harmonics: 



? (EP.) = p^/d 3 X e'El,-x = _i_| d 3 X n^-(iio) 

1 f N 

j^yj d3x II E ^"47ri Ln (p n x)y; nifen (p n )y L 



(X) 



TL — 1 L n .kn 



where X is the unit vector of X. Plugging eq. (110) into (98), we are left with 
an integration of the product of three spherical harmonics having the same 
versor p n as argument. Switching to polar coordinates: 

d 3 p„ — ► dp n p^dfip n and d 3 X — ► dX X 2 d% 

we have [12]: 



(111) 



<2L n + 1 

47T 



(0,0|/ n L n | l n , 0)(jn n , k n \l n L n \l n , Tn n ) 



Furthermore, by recalling the definition of spherical harmonic: 

N ~ 



iV N In T _|_ -I 

n *wx) = n \l^r*>boMoy 



(112) 



n=l 



n=l 



where 6 and f are, respectively, the polar and azimuthal angles of the unit 
vector X, we can solve the angular integration in f^: 
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an ± n 



121,1 :1 D L avAoy 



71=1 

1 r 2w 



An 



N 

n 

n=l 



/2L n + 1 

47T 



TV 

J n=l 



^ /or i i 

= II V 47r (^i, 0; • • • ; L N , 0|0, 0)(0, 0|Li, fci; . . . ; L N , k N ) (113) 

n=l 



47T 



where we have used the unitarity of R and the last equality is a special case 
of eq. (106). 

Finally, the microcanonical channel weight reads: 
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{Nj} 



dX X 2 



N 

II / d P« Pr 

n=l J 



N 



n=l 



X 



II — ^ — huiVrX)- J v dx n x T 2 t jf n (p n x n )(0, 0|/ n L n |/„, 0) 



Jn^n,^ r n ,L n ,kn I 71—1 



x (J 1; Hi, . . . ; Jat, /xjv| J, A) ( J, A| Ji, ^5 . . . ; Jtv, /Xjv) 
x (Lx, 0; ... ; Ljv, 0|0, 0)(0, 0|L l5 fc l5 . . . ; L N , k N ) . 



(114) 



In the above expression, only N + 1 integrations are left instead of 37V, but 
many series appear which can make the numerical calculation as cumbersome 
as in the plane wave expansion (80). 



8 Classical limit 



We now show that the microcanonical channel weight we have calculated in 
the quantum case reduces to what is expected for a fully classical system. The 
classical limit of eq. (98) (we will not consider the quantum statistics case) can 
be obtained by reintroducing H and then taking the limit fi — > 0. Let us then 
rewrite eq. (98) this way, with the invariant SU(2) measure in the axis-angle 
parametrization: 
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n {Nj} = (2J+l)^ / dQ fi / o 47r d^sin 2 ^3[ A (Rri) 



K 1 



X 



ix„/fi-(p„ — R. (V')Pn) 



A? 



s 4 [P -T.pt. 



n=l 



(115) 



Setting: 



the eq. (115) can be rewritten as: 



ft 



Q {Nj} = (2J+l)g^ / d{l h J** dip sm^D^R^)' 1 ) 

tlD n (R A (ljj)) f f iy n .(p n -R. 1 (^)p n ) 



TV 
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,n=l 



(2tt)* 



/ d 3 p„ / d 3 y n e* 



N 



5 4 P-TPn • 



71=1 



(116) 



In the limit ft — > 0, the integration domain V/H 3 of the Fourier integrals in 
y becomes very large and the integrals tend to a Dirac's delta distribution. 
Hence (p n — > R ii 1 (V')pn) for each p n and this means that the rotation R^ 1 
tends to the identity, or, equivalently ip — > 0,2n, An. Indeed, it can be shown 
that for increasingly small ft, the integrand in eq. (116) develops 4 symmetric 
narrow gaussian peaks in ip with maxima at i/j = e, 2ir — e, 2ir + e, 4tt — e 
with e — > as ft — > 0. Hence, we can reduce the integration on the angle ip 
to the interval [0, 7r] multiplying (116) by 4 and in this interval approximate 
sin-0 ~ *0 and cos^ ~ 1. If v is a generic vector, the vector Rh(ip)v reads: 



Rn(VO v — vcos^ + (n x v) sin^ + (1 — cosVO v ■ fin 



therefore, for small ip: 
thus: 



Rfi(V') — v + ^(n x v) 



Pn - Rft 1 (V')Pn - ^(n X p) 

By using (119) and sin^ ~ ip in (116) we obtain: 



(117) 

(118) 
(119) 



n W>£o (2J+1) 



8tt 2 



7=1 iV J' 



X 



II 7^ trD s -(RM) I d 3 Pn / d 3 y n e^^] 5 4 (p-J2 p n ) 

,n=l \ Z7r ) J Jv / h J V n=l J 



(120) 



29 



Now ti\D 5n (Rn) can be replaced with (2S n + 1) if ip — > 0. On the other hand, 
if J/H and A/ft are macroscopic integer numbers, we can also take the large 
spin limit of the rotation matrix and replace D^R^tp) -1 ) with exp[itph- 3/h] 
where J = Ak is defined as the total angular momentum vector. Therefore: 



^(^-EPnJ- (121) 

We can introduce a classical orbital momentum L by noting that: 

y nj ■ (n x p n .) = -n • (y n . x p nj ) = -n • L^/ft (122) 

Also, one can introduce a the three- vector = and rewriting the eq. (121) 
with d 3 = dQ^dipip 2 : 



(2tt) 3 7 Jv/& 



Q {Nj} ~(2J + l)^ jd 3 e^ 



x- 



(27f) 3 Jv/ka 



d 3 y„ e -^' L "/ ft 



g (25,- + 1)"' 

3=1 
N 



r N r 

n / d3 p« 

.n=l J 



(123) 



where the integration in <fi is performed over the whole space as the contribu- 
tion of the ip interval [n, +oo) is negligible for ft — > 0. Now, by restoring the 
variable x = fty and rescaling <fi — > <f>/h we finally obtain: 



{Nj} 
classical 



(2J + 1) 



7T 



TT (2g 3 - + 1)^ 

Mr ^ 



X 



n p^/MH * 4 ( p -|>) * 3 ( J -J>») ■ (124) 



which is the classical expression of the microcanonical channel weight with 
angular momentum conservation [5]. As expected, in the classical case, this 
conservation implies a simple 5 3 as though the angular momentum components 
were commuting observables. To be noted that the contribution of spin reduces 
to an overall degeneracy factor. 
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9 Grand-canonical limit 



It is also interesting to derive the (grand-) canonical partition function with 
angular momentum conservation as a limiting case of the microcanonical par- 
tition function for large volumes and energies. Again, we will confine ourselves 
to the case of Boltzmann statistics. One can start from the expression (98) 
of the microcanonical channel weight and disregard the factor (0|Py|0) as 
P v - I: 



I 

3=1 iy r 



x 



n {Nj} = (2 j + 1) JL J d n h £ ^ sin 2 1 ^(r. (^)-i) 



(125) 

where we have used the axis-angle parametrization. We have seen in the pre- 
vious section, dealing with the classical limit, that for large volumes, the in- 
tegrand in ip in the above equation develops 4 symmetric gaussian peaks in ip 
with maxima at ip = e, 2n — e, 2n + e, 4n — e with e — > as V" — > oo and that 
the Fourier integrals be approximated by: 



Fi/fn -R-(V;l _1 n ) ~ - f d 3 x e ix " i/ ' (fixpn) 

^MPn K nW Pn j ( ^ )3 y y a x e 



(126) 



In view of eqs. (117), (118). Thus, plugging (126), into eq. (125) and using the 
same approximations on the group integration as in Sect. 8: 



j=l 3 



X 



fl / d 3 p„7^/ d 3 x n e-^^ x P") trD s -(RM)} 5 4 (p-^Pr) 

_n=l Jv ( 27r ) JV J V n=l J 



(127) 



A quick way of deriving the grand-canonical limit of the expression (127) is to 
replace the 5 factor of four momentum conservation with the canonical weight 
l\n=i e ~ £n ^ T = Tlf=i Hnj=i e ~ £ " J ^ where T is the temperature. In fact, this 
procedure corresponds to replace the projector 5 4 (P — P), with the operator 
exp(—H/T) so as to match the definition of grand-canonical partition function 
(with vanishing chemical potentials) with fixed angular momentum: 

Zj = tr[exp(-H/T)P JX P v ] (128) 
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Indeed, it can be seen from eqs. (15), (21) that (128) is obtained from the 
microcanonical partition function definition replacing the projector onto fixed 
energy-momentum with the familiar exp[— H/T]. We could have also devel- 
oped an integral expression of the microcanonical partition function summing 
up all channels (127) and inserting a Fourier expansion of the 5 A in order to 
calculate the grand-canonical limit through a saddle-point expansion [6], but 
this would imply lengthy calculations. So, with the above replacement, we 
obtain a grand-canonical channel weight: 



'{Nj} 
K 



X 



(2J+1)^ J dSh£d1; ^^(Rfity)- 1 ) 
II II / / d 3 x % e-^K-^) trZ^'(R^)) e~^ T 



= (2J+1)— IdShi d^^U^R^)" 1 ) 



x 



K -I 1 ^ 

ft 7A3 / d3x / d ' 3 P e-^/ T e-^( xxp ) trL^R^)) 

7=1 \^ ) V 



(129) 



It is now straightforward to calculate the full grand-canonical partition func- 
tion by summing over all possible channels, i.e. over all multiplicities {Nj} 
from to oo. The grand-canonical partition function Zj with fixed angular 
momentum turns out to be an SU(2) integral of an exponential function: 



Z J= 12 ■■■ 12 Z {Nj} 

N 1= N K =0 



= (2j + i)-L f an h ^ 2 D{ x {R h {^)- 1 ) 

Q7T Z J JO 



x exp 



EtA^ / d 3 x / d3 P e- £ ^ T e-^ A '( xx P) tiD s * (R A (iJ>)) 

{27T) 6 JV J 



(130) 



Defining the vector <f> = tph (so that cf) = n, the eq. (130) can be rewritten as: 



Zj={2J+l)^-J d 3 ^(R^)- 1 ) 

J2 7A3 / d 3 x / d 3 P e- £ ^ T e-^( xxp ) trZ^(R^(0)) 

j' = l V 



x exp 



(131) 



For large values of J and A, we can replace the matrix element in (131) with 
exp[i(f> ■ J] where J = Ak, like in Sect. 8 and write Zj as: 
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Zj = (2J+l)-^ [ d 3 e^ J 
8n 2 J\4,\<tt 

K 



x exp 



it-pr^ [ d3x / d3 P e- £ ^ T e-^( xxp ) trZ^(R^(0)) 



(132) 



which corresponds to a classical limit were not for the presence of the traces 
of matrices D s ? '(R^(0)). 



10 Microcanonical ensemble with fixed parity 



In this section we will show how to enforce a fixed parity for an ideal relativistic 
quantum gas. 

As has been shown in Sect. 2, the projector P; onto irreducible Poincare states 
can be factorized as in eq. (14): 

P, = Pp«^ (133) 

where II is the parity of the system. Rewriting the general state weight with 
the full projector: 

"/ = (flPpJX-^-Pvlf) = \(f\PpJxPv\f) + ^(f\PpjxflPv\f) (134) 

we are left with the calculation of two terms: the first is simply eq. (70) times 
a factor 1/2, while the second is a new one to be calculated. To do that, one 
should find the action of 11 on a single particle state. Indeed, there is no unique 
answer to this question, as any definition of the unitary representation of space 
inversion on Hilbert space fulfilling commutation rules of 10(1,3)^ is equally 
valid. This leaves some freedom and we choose: 

f\\p,a) = El n P»^W n P]" ln M) • (135) 

where lip = (e, — p). To be a good one, the definition (135) must meet the 
following requirements: 

• D 5 ([rip] _1 n[p]) must be unitary; 

• for each state |p°,<r), where p° = (to, 0), the space inversion ll must yield 
n|p°, a) = \p°, a)r) where 77 is the intrinsic parity. 

With regard to the first point, it should be noted that D is a finite dimensional 
representation of the universal covering group of SO(l,3)+, i.e. SL(2,C) and a 
unitary representation of the subgroups SO(3) and SU(2) respectively. For it 
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to be extended so as to include space inversion operator, one should recall the 
second Shur's lemma: 

In a finite- dimensional irreducible representation of a group G, the only 
elements which commute with all others are multiples of the identity. 

Since space inversion commutes with all rotations, i.e. [11, R] = 0, then: 

D S (U)D S (R) = D S {UR) = D S (RU) = D S (R)D S (U) 

thus D S (U) commutes with all elements of a finite-dimensional representation 
of a group and, according to the second Schur's lemma, must be a multiple of 
the identity: 

D s (n)=r]l (136) 

Moreover, since D s (r\ 2 ) = D s (\) = if I, then rf = 1 =>• rj = ±1. In the 
second step of the proof, we show that [n£>] _1 n[p] does not involve any boost, 
hence it belongs to the subgroup 0(3). This can be done by noting that, when 
acting upon the unit time vector t = (1,0): 

[np]- 1 ^ — t =► [npj-^tp] g 0(3) . (137) 

Since the above transformation belongs to 0(3), it is either a rotation or the 
product of a rotation by 11 itself. In both cases, because of (136) and rf = 1, 
the matrix .D 5 ([np] _1 n[p]) is unitary and the first requirement is met. 

As far as the second requirement is concerned, noting that [p°] = [lip ] = I 
and using (135), (136): 

fibV) = El n P°^>^([ n /]" ln [p ]) = Eb°.^( n ) 
=Eb°^)^ = b°,<7)j7 ( 138 ) 

where r\ is the intrinsic parity. Hence, the second requirement is also met. 
We now have to check the consistency of (135) with commutation rules: 



nR|p,<j) = RU\p,a) 

f\i\p,a) = L- X n|j9, ct> V|p,cr). (139) 

f\P\p,a) = -Pf\\p,a) 

R being a rotation, L a pure Lorentz boost and P being the momentum op- 
erator. The last equality is trivially fulfilled. For the first commutation rule, 
going back to definition (135): 
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nk\p,a) = Y / n\Rp,T)D s Ta (lRp]- 1 R\p}) 

T 

= ^|nRp, /U ) J D^([nRp]- 1 n[Rp]) J Df (7 ([Rp]- 1 R[p]) 

= ^|nRp,/i) J D^([nRp]- 1 nRb]) (mo) 

This relation tells us that the composition of operators transforming a state 
with a Wigner- rotation's rule (like space inversion itself according to (135)) 
yields an operator fulfilling the same rule. Hence we have: 

Rf\\p,a) = J2\RUp, f ,}D^([Rnp]- 1 RU[p]) (141) 

and taking into account that [l~l, R] = it follows from eqs. (140) and (141) 
that riR|p, cr> = Rfl|p,<7) V|p,cr). 

Similarly, for Lorentz boosts: 

fi i\p, a) = E 1 n Lp, fj) d%{ [n ip] 1 n i[ P ] ) 

L-^b.a) = ^|L- 1 np, /U ) J D^([L- 1 n P ]- 1 L- 1 n[p]) (142) 

thus, taking into account that ULp = L _1 rip also the last equality in (139) is 
proved. 

The conclusion is that the (135) is consistent with all properties needed by a 
good definition of space inversion representation. 

We can now turn to calculate the second term in eq. (134). We will consider 
the special case of a spherical system, so that [11, P v ] = 0, and considering the 
simple case of a single particle state: 

(p,cr\P PJX nP v \p,a) (143) 

Since [11, P v ] = 0, we can rewrite (143) as: 

E (P^ I P ^a Pv I Up, r ) D s Ta ( [ Up] - 1 n [p] ) (144) 

T 

Inserting a resolution of the identity: 

E / dV(p^|PpjAb , ^)(p , ^|Py|np,r) J Df (7 ([np]- 1 nb]) . (145) 
and, by using (28) and (49) we can write: 
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(p,a\P PJX nP v \p,a)=5 A (P-p) |d 3 p' (2J+1) / dR D A J A (R-> 3 (Rp' - p) 

xE^(N" lR b1)J (^(M _1 [np]) + ^([p / ] t [np] f - 1 )) ^([np]- 1 ^]) 

xF y (p + p')(0|P y |0> (146) 
The matrix products in (146) can be worked out as follows: 



E^([ R pr lR M)^([prMnp])^([np]- 1 n[p]) = J Df (T ([p]- 1 Rn[p]) 

= E^([RpT^ 1 b'] t - 1 )^([p / ] t [np] t - 1 )^([np]- 1 n[ P ])t- 1 
= E^(['V] tR [ n p] t " 1 ) (^ 5 ([np] t ) J D 5 (n)t- 1 J D^([p]t-i)) ra 

T 

= E^([fV] tR ) (^ 5 (n)i> fl (|p] t - 1 )) w 

T 

= Di([p]tRn[p]t- 1 ) (147) 

where the constraint Rp' = p has been used as well as the unitarity of the 
Wigner rotation and of the matrices -D 5 ([np] _1 n[p]) and D S (U). Therefore, 
taking (147) into account, the eq. (146) can be simplified into: 

(p,a\P PJX f\P v \p,a)=5 A (P-p) jd 3 p' (2J+1) J dR ^(R-^Rp' - p) 
xi(D 5 (b]- 1 Rn[p])+ J D 5 ([p]tRn[p]t- 1 ) (TCT F y (p + p')(0|Pv,|0) (148) 

If we now sum over polarization states a, the last factor in previous equation 
yields the trace of _D 5 (Rn). Since: 



tr£> 5 (Rn) = tr [D s (R)D s {n) 
we finally get: 



= tr 



D s (R)r]I =r]trD s (R) . 



(149) 



E(^|PpjAnPvb,<x> = ^ 4 (P -p)(2J +1) / dR d^r- 1 ) 

cr J 

x F y (p + R- 1 p)trD 5 (R)(0|Py|0) . 



(150) 



This is, up to a factor r], exactly the same result one would obtain without 
space inversion, the only difference being in the argument Fourier integral 
which is now (p + R _1 p) instead of (p — R _1 p). Thus, it is not difficult to 
realize that the contribution to the microcanonical weight of a generic channel 
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with N different particles, of the space-reflected term in the right-hand side 
of (134) coincides with eq. (61) up to a factor II Y\^ =1 r)j 3 = IHlj^} and a 
"plus" sign in the argument of Fourier integrals. Finally, we quote the full 
expression of the microcanonical channel weight for a spherical system at rest 
in the Boltzmann statistics: 



1 f 47T ib ( 1\ ^ r 

{Nj} = ~J Q ^ sin - sin { J + - j V I] . ^ II J d 3 Pn, 

x \ (F v (p nj - RaW'Pn,) + nn {Jv . } F v (p n . + R3(^) _1 Pn,)) 



X 



sin(5 , i + 



sin(f) 



N 



\Pv\0) 



(151) 



n=l 



which can be obtained by removing all permutations but identity in eq. (96) 
and adding the parity-transformed term. The microcanonical partition func- 
tion can be obtained by summing (151) over all possible channels. 



11 Conclusions 



We have studied the microcanonical ensemble of a multi-species ideal rela- 
tivistic quantum gas fixing the maximal set of observables pertaining to the 
space-time symmetry group. In the rest frame of the system, where P = this 
means enforcing the conservation of energy-momentum, total angular momen- 
tum and parity. We have provided a full consistent treatment and solved the 
problem for the most general case of particles with spin, generalizing the re- 
sults of Cerulus [2]. 

The microcanonical partition function Q has been expressed as a sum over 
channels, i.e. at fixed multiplicities of each particle species, each term be- 
ing defined microcanonical channel weight. These have been obtained in a 
quantum field theoretical framework, applying to spinorial fields the formulae 
derived for scalar fields in our previous work [1]. This extension relies on two 
reasonable, yet unproved assumptions: 

(1) the validity of eqs. (47) and (48) which is conjectured on the basis of the 
V — > oo limit and the fact that they have been proved for the scalar field 

[i]; 

(2) that also for particles with spin, the expressions of the channel weights 
involving identical particles are the same as those with different particles 
with (anti)symmetrization, that has been proved for the scalar field [1]. 

The implementation of angular momentum conservation leads to an expression 
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of the microcanonical weight of a channel involving a further integration over 
the parameters of the SU(2) group, with respect to the case of only energy- 
momentum conservation. This expression nicely agrees with that obtained 
with an expansion onto the angular momentum basis and leads to the expected 
classical limit, already known and used in literature. We have been working 
out the grand-canonical partition function with fixed angular momentum in 
both the case of small and large angular momentum. We also obtained the 
microcanonical partition function of a system with fixed parity. 
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